Spin-rotation coupling in non-exponential decay of hydrogenlike heavy ions 
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We discuss a model in which a recently reported modulation in the decay of the hydrogenlike 



ions u p r ats + anc j Pm arises from the coupling of rotation to the spin of electron and nuclei 
(Thomas precession) . A similar model describes the electron modulation in muon g — 2 experiments 
correctly. Agreement with the GSI experimental results is obtained for the current QED-values of 



the bound electron g- factors, g( Pr' 



1.872 and g( 14 Tm 



60+ \ 



1.864, if the Lorentz factor 
of the bound electron is ~ 1.88. The latter is fixed by either of the two sets of experimental data. 
The model predicts that the modulation is not observable if the motion of the ions is linear, or if 
the ions are stopped in a target. 

PACS numbers: 23.40.-s, 27.60.+J 



Recent experiments carried out at the storage ring ESR 
of GSI in Darmstadt [H, HJ reveal an oscillation in the or- 
bital electron capture and subsequent decay of hydrogen- 
like 140 Pr 58 + and 142 Pm 60+ . The modulation has a pe- 
riod of 7.069(8) s and 7.10(22) s respectively in the labo- 
ratory frame and is superimposed on the expected expo- 
nential decay. The " zero hypothesis" of a pure exponen- 
tial decay has been excluded at 99% C.L. for both ions. 
Equally improbable causes seem to be periodic instabil- 
ities in the storage ring and detection apparatus. The 
effect has been tentatively attributed to neutrino flavor 
mixing, but this interpretation presents difficulties that 
have been discussed in the literature 0]. 

It may be useful to analyze the effect of spin-rotation 
coupling on the ESR modulation effect. It is in fact 
the spin part of the Thomas precession, or spin-rotation 
coupling, that is responsible for a similar modulation in 
g — 2 experiments. The corresponding modulation in the 
detection of electrons produced in the in-flight decay of 
muons in storage rings yields the anomalous part of the 
muon magnetic moment directly [4j. 

An important point raised in work on 140 p r 58 + and 
i42p m 60+ j g fogfc m t ne i n iti a i state these ions, with nu- 
clear spin 1=1, and the bound electron can have total 
angular momentum 1/2 when the electrons and nucleus 
spins are antiparallel, or 3/2 when the spins are parallel. 
The final state, however, can only have spin 1/2 because 
decay from the spin 3/2 state is forbidden by the conser- 
vation of the F quantum number. 

For the sake of simplicity, we treat the nuclei as spin- 
1 particles interacting with their single electron by an 
amount sufficient to drag it along the same cyclotron 
orbits in the storage ring. We use units h = c = 1, but 
re-introduce standard units when discussing the results. 

Following Bell and Leinaas [j| , we write the full Hamil- 
tonian that describes the behavior of nucleus and bound 



electron in the external induction field B of the ring as 
H = Hq + i? (e) + H {n) , 

where 



(e) , rr(n) 



Ho = + H, 



' 



H (e) = -2g e ^ B s ■ B - u e ■ (L + s) 
HW = -g nl x N l ■ B - w„ • (L + 1) . 



(1) 
(2) 



i/p 6 '™-* contain all the usual standard terms (Coulomb 
potential, spin-orbit coupling, spin-spin coupling, etc.), 
g e and g n are the g- factors, and [Ib and /ijv are the Bohr 
and nuclear magnetons. 

The Hamiltonians and contain the Pauli 

spin matrices via s = <r/2 and the 3x3 matrices I and 
the last terms account for the Thomas precession. If the 
orbits in the ring are stable, the coupling of w to the or- 
bital angular momentum L may be neglected ■ We will 
also neglect, for simplicity, any stray electric fields and 
electric fields needed to stabilize the orbits and all quan- 
tities in {1} and ([2]) will be taken to be time-independent. 

One can also arrive at (p} and ([2]) starting from the 
Dirac equation [fj, @, 0, S] and the corresponding rela- 
tivistic wave equation for spin-1 particles Q- The im- 
portance of spin-rotation coupling has also been greatly 
emphasized, in a different context, by Mashhoon [9(. 

and contribute to the overall energy 

and of the states. By referring (JlJ and |2]) to a 

left-handed tern of axes rotating about the X2-axis in 
the clockwise direction of the ions and with the £3-axis 
tangent to the ion orbit in the direction of its momentum 
, we find that the operators governing the evolution of 
the spin states are 



E' : 



M e = 



fi e B 



E' : 



(3) 



2 



and 



M " 



V 



r 

E n - i- 
2 

-i (U! n - H n B) 





r 

J5 n - i- 
2 





/ 



(4) 

where T is the energy width of the nucleus. Instead of 
P|, we consider the sub- matrix 



/ 



M n 



r 

B n - i- 
2 



i (U n - UnB) 

r I . (5) 



because only the non-zero, off-diagonal terms are rele- 
vant. 

Before the ion decays, the nucleus and the electron 
separately can be represented by [13] 

\i, l {t)) = a l {t)\^ l + ) + b\t)\^_), (6) 

where I = e,n and and are the two helicity 

states. The coefficients a (t) and b l (t) evolve in time ac- 
cording to the equation 



a l (t) 
b l (t) 



(7) 



which can be solved by diagonalizing the matrices M l . 
The spin flip probabilities for electron and nucleus are 



PW = |(^|r(*)>| 2 = ~[l-cosO e t] 



and 



where 



p(n) 



-rt 



i fl n i] 



(8) 



(9) 



fh 



gefJ^B - u e , fi n = 2(ft„B - u. n 



The probability for electron and nucleus to flip their spins 
at time t, if at t = their spins were both antiparallel to 
the nucleus momentum, is 



P oc p( e )p(") = 



,-rt 



I — cos fLi — cos fLt- 



(10) 



cos(i7 e + fi n )t + cos(fi e — fin 



)<]}. 



We use the uncorrelated probability ([TO]) because in our 
model the electron and nucleus spins are treated as inde- 
pendent. The second and third terms in (fT0|) represent 
the spin-flip angular frequencies of nucleus and electron. 
The latter frequency, when applied to the muon, gives 
the muon anomaly in g — 2 experiments. The fourth 



term represents transitions to the forbidden spin-3/2 fi- 
nal state with angular frequency fi' = f2 e + fi n . The last 
term yields the probability that the electron and nucleus 
spins are antiparallel, which is the probability of interest. 
We obtain 

„-rt r 



P' 



1 



1 H — cos(O e — Q n )t 



(11) 



The angular frequency fi = f2 e — fi n with which the elec- 
tron and ion spins find themselves antiparallel and which, 
therefore, enables the ion decay, must now be related to 
laboratory system quantities. 
We now follow 



111 . 1 1 21 ] and re-introduce standard 



units. 

The calculation of the precession frequency of the elec- 
tron spin in its motion about the heavy nucleus first takes 
into account the contribution given by the first term in 
([1]) which in the nucleus rest frame is given by ui g — 

— - — B. We also assume that (f3 e ■ B e ) = = (/3 e x E e ), 
where B e and E e are the magnetic and electric fields in 
the nucleus rest frame and f3 e = v e /c is the velocity of 
the electron relative to the nucleus. The spin-rotation 
term in (TT]) is given, in the rest frame of the electron, by 



WT ,rf 



7e 



7e + l 



e x f3 e , where e is the electron acceler- 



ation, and 7 e = 1 / \J 1 — (3% is the Lorentz factor. o?t ,r/j 
hence e , are generated by the Coulomb force Fc = eE n , 
where E„ is the electric field of the charged nucleus, and 
by the Lorentz force = e(E' + /3 e x B'), where E' and 
B' are the electric and magnetic fields in the rest frame of 
the electron. They are related to the magnetic field of the 
storage ring by means of a Lorentz transformation. Since 
the electron propagates in an external magnetic field, its 
velocity is shifted by ~ [iBBlipiool 2 . This term, however, 
turns out to be negligible. 

The electric field E„ gives rise to the anomalous Zee- 
man effect. This term is not relevant to our study and 
will not be considered further. 

The force Fl, on the other hand, generates an ac- 
celeration of the bound electron that contributes to the 
Thomas precession, hence to the precession of the elec- 
tron spin with respect to the magnetic field B. Pass- 
ing from the electron rest frame to the nucleus frame 
and using the Lorentz transformation of the electric and 
magnetic fields, the Thomas precession becomes ijJt = 
e 2( 7e -l). 



7e 



-B. Therefore, the precession frequency of 
ds , 



the electron is given by ~j~lnf — x s , where 



Q e = U)g + U>T = — 



2ms 



a e - H ) B 

7e 



~ -1.7587 x 10 11 f a e - 1 + — ) —Hz u 2 .(12) 

a e = (|<?e| ~ 2)/2 is the electron magnetic moment 
anomaly, and u 2 is a unitary vector direct along the xi- 



3 



axis. It is parallel to the magnetic field of the storage 
ring and therefore orthogonal to the circular orbit of the 
nucleus. 

If the nucleus moves along a circular orbit in the pres- 
ence of only the external magnetic field, then the preces- 
sion frequency of the nucleus spin, in the reference frame 

comoving with the nucleus, is ^~lnf = x ^ 
where 



II, m 



^-lU^B 

HQ J A m p 



7 B 

~ 4.785 x 10 7 — Hz u 2 . 



2Z J Am p 

(13) 



A and Z are the mass and charge numbers, and we have 
written ji n = fl(eh/2m p ). The spin of the nucleus does 
not precess if p n M/HQ = 1. In (Tl3|) we have approx- 
imated M by Am p . The calculation of M can be im- 
proved by including contributions from nucleon binding 
energies, nuclear surface and Coulomb repulsion effects 
and the Pauli exclusion principle. Not all these contri- 
butions are positive. 

We must now calculate the relative spin precession that 
the motions of electron and nucleus generate. From (fTZ)) . 
([TBI we obtain 



n 



2 

Ye 



2Z 
~A 



eB 



(14) 



The magnetic field is related to the the angular velocity 
of the particles in the laboratory frame 2irf by 

QB _ 7 2 (2W) 2 p 
M Pc 



where 7 = — (3 2 is the Lorentz factor, and p is the 

radius of the orbits. The revolution frequency / of the 
primary beam can be obtained from / = v/L, where L 
is the length of the closed orbit, and from the definition 
of magnetic rigidity Bp = Mvj/Q 2]. We get 



/ = 



Z Bp 



1 



A L 




2MHz. 



(15) 



It is the frequency Q/2n that must be compared with 
the experimental signal OAAHz found for 140 p r 58 + and 
I42p m 60+ jjy means f t ne equation 



fl = 2tt0.14Hz. 



(16) 



The calculation of g-factors, based on bound state (BS) 
QED, can be carried out with accuracy even though, in 
our case, the expansion parameter is Za ~ 0.4. The re- 



sults agree with available direct measurements 13j-[l 



In particular, the BS-QED calculation given in [2(| in- 
cludes radiative corrections of order a/ir and exact bind- 
ing corrections. It yields 



1 + 2^/1 - (azy 



+ -C^(aZ) 



(17) 



where can be approximated by 

C( 2 >(aZ)^ + l( a Z) 2 + I(aZ) 4 . 

The values of g\ calculated by applying (TP?) to C 5+ and 
7+ are in good agreement with the experimental re- 
sults. The same formula (HTl) extrapolated to Z = 59 
and Z = 61, gives g\ ~ 1.87205 for 14 °p r 58 + and 
g b e ~ 1.86365 for 142 Pm 60 +. The addition of more ex- 
pansion terms following the formulae given in flit ] does 
not change these results appreciably. The correspond- 
ing anomalous parts become respectively a e ~ —0.06397 
and a e ~ —0.06817. The only free parameter in (fM]) 
is the electron Lorentz factor j e . It can be fixed by 
using one set of experimental data and using the re- 
sult in the second set. Alternatively, we can choose as 
free parameter the average distance R of the electron 
from the nucleus and estimate the electron velocity by 

, 1 Ze 2 
means ol the equation ra e ^ e v e = 



ft 



b 2 — b ] , where b 



47T£ R 

1 Ze 2 



We obtain 



From 



47T£o m e c 2 R 

the experimental values L = 108.3 m, Bp = 6.44 Tm, 
= 1.43, (3 = 0.71, assuming that Jx — 2.5 as given in 
1|, using the values of a e just found and treating R as a 
parameter, we find that equation (|16p has the solutions 
7 e (Pr 58+ ) - 1.88 1 35 and 7 e (Pm 60+ ) - 1.87392, which 
correspond to R ~ 123 fm ~ 0.145 az for Pr 58+ , and 
R ~ 128 fm - 0.151 as for Pm 60 +. Here a z {= a /Z) is 
the Bohr radius (ao ~ 0.53 x 10~ 10 m). As a comparison, 
the nuclear radius is given by Rn ~ R0A 1 / 3 fm ~ 6.2 fm, 
in nuclear models in which the nucleons are uniformly 
distributed in the nucleus (Rq ~ 1.2 fm). 

In our calculations we have adhered to the assignment 
p, = 2.5 [l|. However, a need for a re-measurement of 
nuclear magnetic moments arises from improved ways to 
probe QED effects, as pointed out in 17 1. 

Eq. ([TJ} does not depend on h. This should be ex- 
pected because the electron orbits and those to which 
(fl~5|) refers are classical. This is also the extent to which 
the treatments of particles in storage rings given in [j| 
and 11 1 agree (see the discussion given in Q). Bell 
and Leinaas consider quantum fluctuations away from 
the classical orbit. We neglect them, in agreement with 
Jackson [111 ]. 

The model has consequences that must be emphasized. 
Identical ions in rings of different radii p\ and P2, but 
with the same magnetic rigidity, and the same speed, give 



4 



rise to Q ~ Therefore, the effect becomes smaller in 
larger rings for fixed A/Z, g e and fi. 

It is also clear that according to our hypothesis, which 
is essentially based on the coupling between rotation and 
the spins of electron and nucleus, the GSI effect disap- 
pears when / = 0. Accordingly, the effect can not be 
measured in experiments based on the linear motion of 
the ions, nor can it be observed by stopping ions in thick 
aluminum foils [HI • 

A more complete model would also require the rela- 
tivistic treatment of hydrogen-like atoms (with large nu- 
clear charge) in a noninertial reference frame. 

In summary, the GSI experiments measure fl in our 
model and are well suited to study nuclear radii, g e and 
fl(g e )- This is important in fields like nuclear physics, 
QED, BS-QED and stellar nucleosynthesis (2fT |. 
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